CHAPITRE 3 



EQUATION DE REYNOLDS 

1 - INTRODUCTION A LA LUBRIFICATION 



La lubrification hydrodynamique est un domaine important de la tribologie, c’est 
l’êtude des contacts dans lesquels un film de fluide sêpare les surfaces en prêsence. Dans le 
cas ou le film de fluide sêpare totalement les surfaces, les aspêritês et les dêfauts de fonne ont 
des dimensions infêrieures a l’êpaisseur du film. 

le problème peut se schêmatiser de la fagon 

Un film de fluide visqueux sêpare les deux 
surfaces du mêcanisme ; il s’agit de calculer : 
la charge W que peut supporter le contact, 
la force F ou le couple de frottement (? , 
le dêbit Q du fluide dans le mêcanisme, 
la puissance V dissipêe dans le contact. 
Ainsi il faut dêterminer le champ de pression 
et le champ de vitesse dans le fluide. Par 
ailleurs, l’intêgration des contraintes de 
cisaillement aux surfaces du contact donnera 
la force ou le couple de frottement. 

La pression peut être crêêe par une pompe extêrieure au contact c’est l’hydrostatique : 
c’est le cas par exemple de certains paliers et des glissières de certaines machines outils. 

La pression peut être crêêe par le dêplacement relatif des surfaces, c’est 
l’hydrodynamique : cela conceme les butêes, les paliers fluides, les garnitures mècaniques. ... 

2 - APPROCHE PHYSIOUE DE LA LUBRIFICATION HYDRODYNAMIOUE 

2.1 - Types d’ecoulement dans le film 



Dans tous les mêcanismes êtudiês, 
suivante (fig. 1) : 




Fig. 1 : Schèma d'un contact lubrifiê 



Considêrons, (fig. 2), l’êcoulement plan entre deux surfaces entièrement sèparèes par 
un fluide visqueux newtonien. Si dans le film lubrifiant, les forces massiques extêrieures et les 
forces d’inertie sont nêgligeables devant les forces de viscositê et de pression, l’êcoulement 
rêsulte de la supei'position de deux êcoulements simples : 

L’êcoulement de Couette du au mouvement des surfaces, (fig. 2a). La variation de vitesse 
est linêaire suivant l’êpaisseur du film. 

L’êcoulement de Poiseuille du a un gradient de pression dp/dx dans le sens de 
l’êcoulement. La variation de vitesse suivant l’êpaisseur du film est parabolique, (fig. 2b). 

On peut ainsi rencontrer diffêrents cas de figure, selon le signe des vitesses Uj , U 2 et 
du gradient de pression dp/dx (fig. 2). 
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a) Ecoulement 
de Couette 



b) Ecoulement c) Ecoulements rêsultants 

de Poiseuille 

Fig. 2 : Diffêrents types d’êcoulement 



2.2 - Explication de la portance dans le coin d’huile 



Soient deux plaques planes non parallèles sêparêes par un film fluide newtonien, 
(fig. 3). La plaque infêrieure est animêe d’un mouvement de translation de vitesse ET 
constante. 

Dans le cas d’un êcoulement unidirectionnel 
(problème plan) et s’il n’existe pas de 
variations de pression dans le film, la 
distribution de vitesse dans le contact est 
linêaire (êcoulement de Couette), fig. 3. Dans 
ces conditions le dèbit, pour une largeur L 
selon oz est : 

A l’entrêe Q e = ^ ^ 




A la sortie Q s = 



2 

Lh, U 



Fig. 3 : Hypothèse ; êcoulement de Couette 



Pour un fluide incompressible le dêbit 
volumique doit être constant. Soit : Q e = Q s . 

Ceci est impossible car h e est diffêrent de 
h s . II y a donc nêcessairement gênêration de 

pression dans le film. Sous l’action du 
gradient de pression, le fluide est ralenti a 
l’entrêe et est accêlêrê a la sortie, (fig. 4) de 
telle sorte qu’il y ait conservation du dêbit. 




Fig. 4 : Ecoulement dans un coin d’huile 



2.3 -Relations elementaires 



Considêrons deux surfaces planes, 
(fig. 5) entre lesquelles s’êcoule un fluide 
newtonien. On notera U la vitesse de surface 
infêrieure. L’êcoulement est supposê 
laminaire, permanent, isotherme et 
unidirectionnel. Etudions l’êquilibre d’un 
parallêlêpipède êlêmentaire de côtês dx et dy 
et de largeur unitaire selon Oz, (fig. 5). 



K- 




h(x) 



I I * 

dx 

u 






Fig. 5 : Equilibre d’un parallêlêpipède 
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Le parallêlêpipède est sollicitê par : 

Les contraintes de pression, si la distance h qui sêpare les 2 plans est suffisamment petite 
et si les variations d’êpaisseur du film selon 0x sont faibles, on peut supposer que la 
pression est constante selon l’êpaisseur du film. 

Les contraintes de cisaillement. 

La condition d’êquilibre s’êcrit : 

[p(x)-p(x + dx)]dy + [r xy (y + dy) -x xy (y)]ix = 0 



avec d’après la loi de Newton x = q— ; ou u reprêsente la vitesse du fluide selon 0x et q 



dy 



sa viscositê supposèe constante. Un dêveloppement limitê au premier ordre permet d’êcrire : 

I dx dy = 0 



dp , , d ( du) 

dx dy + — u — 

dx d y 3 y 



soit 



dp 

dx 



[X 



a 2 u 

dy 2 



Cette èquation est une forme très simplifiêe de l’êquation de Navier Stokes en 
êcoulement plan. Compte tenu des hypothèses effectuêes prêcèdemment la vitesse s’ècrit : 



j_dp , 

2\x dx 



+ Cj y + C 2 



ou Q et C 2 sont deux constantes dêterminêes par les conditions aux limites sur la vitesse. 
S’il n’y a pas glissement a la paroi on a : 

pour y = 0 , u = U 

pour y = h , u = 0 



II vient : 



1 

u = — 
2 q 



dp 

dx 



y (y - h) + 



h-y 

h 



U 



Dans une section droite quelconque le dêbit q(x) par unitê de largeur est donnê par la 
relation : 



q(x) = 



Q(x) 

L 



,h(x) h dp hU 

u dy = + 

0 12 p, dx 2 



Pour un fluide incompressible et un êcoulement unidirectionnel le dêbit q (x) doit être 
constant. Soit : 
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h 3 dp hU 

+ 

12p dx 2 



= Cste = 



h U 
2 



ou h est l’êpaisseur du fdm pour laquelle le gradient de pression s’annule. Soit encore : 



dp 

dx 



= 6pU 



h-h* 

h 3 



Cette êquation est une forme très simplifiêe de l’êquation de Reynolds ; elle permet 
d’êtudier des êcoulements plans dans le cadre d’hypothèses simples. 

3 - EOUATION DE REYNOLDS 

3.1 - Derivation de l'equation de Reynolds 

Les êquations gênêrales de la lubrification hydrodynamique par fluide newtonien 
prennent diffêrentes formes selon les hypothèses effectuêes. Dans ce chapitre, nous 
dêmontrons une forme simplifiêe de ces êquations, appelêe êquations de Reynolds [1 a 5]. 

L’êquation de Reynolds est dêduite des êquations de la mêcanique des milieux 
continus, appliquêes a un fluide newtonien et qui peuvent s’êcrire sous la forme suivante : 

la loi de conservation de la masse : 



d p 
ôt 



<3 

+ T (P u i) = 0 

O X; 



la loi fondamentale de la dynamique : 



/ 

P 

V 



du; 
d t 



+ U J 






d x 



J / 



= P f i + 



d°'j 

ÔX: 



la loi de comportement rhêologique, qui pour un fluide newtonien s’êcrit : 

a ij ■= (-p + ^-0)ôij + 2 ps - 



Dans ces expressions, X;et t sont les variables d’espace et de temps, u ; les 
composantes de la vitesse, p la masse volumique du fluide, p la pression, f ; les forces 
massiques extêrieures, 0 - le tenseur des contraintes donnê par la loi de comportement 
rhêologique du fluide, ò- le symbole de Kronecker ( òy = 1 si i = j ; òy = 0 si i ^ j) et ou 
L et p , coefficients de viscositê de Navier, sont des fonctions de l’espace et du temps. 

, 1 , . ,, . ^ du, du 7 du, 

0 , le taux de ddatation cubique s ecrit : 0 = + + 

ôXj dx 2 dx 3 



£jj , le tenseur des taux de dêformations est donnê par : £ ;j 



2 



' du ; 
ôx ( 

V J 



+ 



dUj N 
ôx ; , 
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En reportant l’êquation de comportement rhêologique dans l’êquation de la dynamique 
et en nêgligeant les forces massiques extêrieures, on obtient les êquations de Navier qui 
s’êcrivent : 



ÔU; du ; 

+ u , 



d t 



d x 



J J 



dp d 2 Uj 

+ A — + q 

dx ; dX;dXj 



/ a 2 



Ô“U; d Uj 

+ 



ÔXj dx { d Xj 



9Uj d'k ( du, ôUj\ 



ôXj dx { 



V ôx j 3x ./ 



d q 
d x ■ 



Dans ces relations, i correspond a la direction considêrêe et j joue le rôle d’indice 
muet. En lubrification, la dimension suivant l’êpaisseur du film (0,x 2 ) es t t res faible par 
rapport aux deux autres, ce qui permet, pour dêterminer l’ordre de grandeur des diffêrents 
termes de l’êquation de Navier, d’êcrire [6] : 



x 



x, = 



1 



L 



5 A 2 



X 0 = 



- Xi - tV 

. x 3 -~T, t = — 



L 2 77 



u 



U, = 



1 7 T 



V 



Ui = 



, U 2 



H L L 

u 2 L u, 

VH 2 V 



- T k 

p = — et K = — 

fio 



Dans ces expressions L et V sont respectivement la dimension et la vitesse 

VH 

caractêristiques selon les axes (0,Xj) et (0,x 3 ), H et — j— - celles suivant (0,x 2 ) ; 

— , q 0 et A 0 dêfinissent l’ordre de grandeur du teinps et des coefficients de viscositê. Nous 
supposerons que ces deux coefficients sont du même ordre de grandeur. Ce changement de 

H 2 

variable conduit a poser p = p . Remarquons que la vitesse caractêristique V peut être 

fioVL 

dêfinie de plusieurs faqons ; en hydrodynamique, V sera la vitesse d’une des parois du 
contact. Si les parois du contact sont immobiles, le contact est hydrostatique, on dêfinit alors 
la pression sans dimension par p = p/P s ou P s est, par exemple, la pression d’alimentation ; 



la vitesse caractêristique peut alors s’êcrire : V = 



P S H 

q 0 L 



Ce changement de variable met en êvidence les 3 nombres sans dimension suivants 

le nombre de Reynolds relatif a l’êpaisseur de l’êcoulement, 
fio 

H _3 

e = — le paramètre d’êchelle, qui en lubrification est de l’ordre de 10 , 

A 0 

— le rapport des coefficients de viscositê, gênêralement de l’ordre de 1 . 
fio 



Compte tenu de ces changements de variable l’êquation s’êcrit : 
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dp 
9 X, 



= 8 



-91 



^du, _ 9u, ^ 

A + u 



\ 



9 t 



J 9 x 



+ 8 



J 7 



(- Y X o\ 
|+ + A 

Po 



9 x. 



( £V 

\ ax i/ 



+ 8 [X 



(d 2 u x 9 2 u^ 

+ 



Ôx 2 



Ô X 3 y 



+ [X 



d 2 u 



Ôx 2 



1 + 8 2 



dk d Uj ^ 2 3jx chiy t 9p / 9ui ^ 9u 3 ^ 



|x 0 9 Xj 3 Xj 



ô X[ 3 X[ 3 x 3 I ô x 3 



d x 



>7 



I dp f , c 2 ÔU 2 ^ 



ô x 2 l 9 x 2 



ô x 



i7 



^_ e2 

9 x 0 



-9t 



A u 



9 t 



2 + u 



ÔU 



TT \ 



J XI T7 



9 X: 



+ 8 |X 



( d 2 U 2 

9 x 2 



a 2 u 2 x 
ôx 2 



/ 



p:+^ 

Po 



9 (dUj) 



dx j 



9 x~ 



\ J / 



_ d u 2 k 0 9 X / $ u j 



|x 



Ô X 2 



Lo 



ôx j 



+ 2 - 



Ôp,ÔU 2 Ôp/dU[ 2 9 U 2 \ ô jx / d u 



9 x 2 ô x 2 



9 x. 



9 x~ 



+ 8 ‘ 



9 x 



17 



9x 2 



9 x- 



-+ 8 ‘ 



9u 



TT \ 



9x, 



9 p 
9x, 



= 8 



-9t 



/ 9 u 



9 t 



3 +u 



9u ^ 



J 9 x ■ 



+ 8 



J 7 



/_ r X 0 \ 

k Ll + A 

Lo 



9 ^ ô xxj ^ 



v« x ,/ 



9 x, 



+ 8 p 



/ 9 2 u 3 9 2 u 3 ^ 



9 xf 9 x 3 



+ p. 



■\ 2 - 

. 9 u 



9 x 2 



3 2 

+ 8 



^o 9X c> u j ^ 2 9 pt 9 u 3 ^ 9[x /9 u 3 ^ 9u t ^ 



p 0 9 x 3 9 Xj 9 x 3 9 x 3 9xj l 9xj 9x 3 



9 p. 



9 x- 



/ 9 u, 
9 x, 



+ 8 ‘ 



9u 



TT \ 



9 x, 



En nêgligeant les termes multipliês soit par e 2 , soit par e 3 9t il vient : 



9 p 

9 X; 



Les termes en e 9t , ou termes d’inertie, de l’êquation ci dessus seront nêgligeables 
devant les termes de viscositê quand e 9t « 1 ; ce qui est souvent vêrifiê dans les 
êcoulements rencontrês en lubrification. En effet, dans le cas du contact d’une bille sur un 
plan par exemple (reprêsentatif du roulement a billes ou de l'engrenage) : 

L = 10“ 3 m, H = 10“ 6 m, V = 10m/s, q 0 = 0,05Pa.s, p 0 = 900kg/m 3 

£ = 1,8 10“ 4 

Dans le cas d’un palier a huile de dimensions moyennes : 

L = 0,05m, H = 10“ 4 m, V = 30m/s, p 0 = 0,05Pa.s, p 0 = 900kg/m 3 

e 91 = 1,08 10“ 2 

Des exceptions se rencontrent toutefois quand la viscositê du fluide est faible, quand 
l’êpaisseur du film est importante ou pour les paliers de grandes dimensions. Ainsi pour le 
palier prêcêdant mais lubrifiê avec de l’eau on aura : 



= -891 



9 U: _ 9u 

r- + U : 

9 1 



J 9 x 



j 

9p 

9 x 0 



+ p 



9 2 u, 

9x 2 



9u ; 9 |x 
9 x 2 9x 2 



(i = 1 et 3) 



= 0 
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m = lmPa.s 



et 

e9i = 0,6 



p 0 = 1000kg/m 3 



De même pour un palier de grandes dimensions tel que ceux supportant les lignes 
d’arbre des centrales êlectriques modernes. 

L = 0,5 m, H=10" 3 m, V=100m/s, q = 0,05 Pas p =900 kg/m 3 

s9i = 3,6 

Dans ces cas les forces d’inertie ne sont plus nêgligeables. Ce problème sera abordê 
ultêrieurement. Si par hypothèse on nêglige les forces d’inertie il vient : 



d p _d 2 U; ôUj d u 

— — = U C + 

d x : 



d XŌ u /V 2 u a 2 



1 " 2 ÔX, ()x- 

dp 



(i = 1 et 3) 



d x- 



= 0 



et en revenant aux variables dimensionnêes en coordonnêes cartêsiennes (0 x y z) : 

Op d I ôu) 
d x dy ^ d y 



^P-o 

;>y 

dp d l dw\ 
dz d y ^ d y 



Pour obtenir le champ de vitesse, il faut intêgrer. Si u est indèpendant de y, il vient : 



soit : 



p 

p 



d U d p 

= — -y + A 

d y d x 

d w d p _ 

= — y + C 

d y d z 



1 d p 

u = 

2p d x 

1 d p 



w = 

2p d z 



2 



y 



2 



y 



+ — y + B 

p 

C 

+ — y + D 

p 



Les conditions aux limites sur la vitesse permettent de dêterminer A, B, C et D. Si l’on 
suppose qu’il n’y a pas de glissement du fluide aux parois, elles s’êcrivent, (fig; 6) : 



fy = 0 u = U ^ ; v = 0 et w = W^ 

1y = h u = U 2 ; v = V 2 et w = W 2 
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Fig. 6 : Système d'axes en coordonnêes cartêsiennes 



On a donc : 



Uj =B 

U = J_^P_ h 2 + Ah + B 



2[i d x 






soit : 



A U, -U, 1 dp 

— h 






h 2 [x 5 x 



il vient alors : 



de même pour w : 



1 dp Uj-Uj 

u = — — y(y-h) + y + U^ 

2 p, d x h 



1 ôp . W 2 - Wl 

w = — — y(y-h) + ^- y + Wl 

2 p, d z h 



En reportant ces vitesses dans l’êquation de continuitê et en intêgrant a travers 
l’êpaisseur du film l’inconnue v disparait. En effet, l’êquation de continuitê s’êcrit : 



dP 3 d d 

— + — (pu) + — (pv) + — (pw) = 0 
dtdx dy a z 



en intêgrant, on obtient : 



/o 77 d y + / 0 h f (P u )dy + / h ^-(pv)dy + / h ^-(pw)dy = 0 
JU d t J d x Ju dy J dz 



II vient ainsi : 



f (p v) dy = [p v J = p V 2 car V, = 0 par hypothèse. 

J ° d v 



dy 



Par ailleurs, si p n'est pas fonction de y, on peut êcrire : 
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d p 
d t 



Pour intêgrer les deux autres termes, il faut intervertir les signes dêrivês et sommes en 
remarquant que : 



h(x,z,t) 3 F(x, y, z, t) 
dx 



, Ô »h(x,z,t) 



F(x,y, z, t) dy - F(x,h, z, t) 



d h 
d x 



Soit 



r h d d T T d h 

f — (pu)dy = — f pudy-pU 2 — 
J0 d x d x J0 d x 

h 3 5 h ô h 

r ~ ( p w ) dy = ~ r p w dy ~ p w 2 — 

J0 d z d z J0 d z 



Ces deux intêgrales se calculent en utilisant les expressions donnant les vitesses dans 
le film. On a ainsi : 



soit encore : 



d 


r h 


_p_ a_p 


dx« 


1 0 


2p dx 


d 


■ 


ph 3 dp 


d x 




12p dx 



U2-U1 

h 



y + p U, 



dy = 



-h + pUj h 



de même : 



d /.h 

J^pucly- 



d X 



d /J 1 

Jpwdy: 



d z 



d 


'ph 3 dp^ 


d 

_l_ 


(p u * +u l 


d X 


^ 12p dx j 


d x 


1 P 2 


d 


Oh 

+3 

Q. 


d 1 

_i_ 


' w, + w, 


d z 


^ 12p d z j 


3 z | 


p 2 

\ 



h 



h 



et en regroupant tous les termes il vient : 



d 


^ph 3 3p) 


d 

_i_ 


^ph 3 3p) 


d | 


( P (U 2 + Uj)h\ 


d x 


^ 12p d x f 


dz 


^ 12p d z ^ 


d x 1 


l 2 J 



d f p(w 2 + w t )h 

dz\ 2 



d h 

- pU 2 — 
d x 

d h , d p 

- pw 2 — + p V 2 +h — 

d x d t 
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soit encore : 



d 


m 

j=i 

Q. 


d p ' 


d 

_i_ 


Vh 3 


d p ' 


d x 


{ 0 


d x^ 


dz 


1 * 


d 2 ] 



d h d h 

= 6p(U, - U 2 ) — + 6p (Wj - W 2 ) — 

d x d z 



+ 6h — — — [p (Uj + U 2 )]+ 6h — [p(W, + W 2 )]+ 12 pV, +12h 
d x d z 



d p 
d t 



A cette êquation, dite êquation de Reynolds, s'ajoutent les êquations donnant les 
vitesses dans le fluide : 



1 d p 

u = 

2p d x 

1 d p 

w = 

2p dz 



y(y-h) + V^u l+ 7:U, 

h h 

.y(y- h ) + ^ Wl +^W 

h h 



2 



Ainsi que les contraintes de cisaillement dans le fluide : 



x 



xy 



X 



y z 



= p 



= p 



du 
d y 
d w 
d y 



V (2y-h) + (U, -U,A 
2 d x h 

j-^(2y-h) + (W 2 -W,)^ 
2 dz h 



En coordonnêes cylindriques, (fig. 7) une dêmarche analogue permet d'obtenir les 
êquations suivantes. Dans ce système d'axes (O r, 0 , z) les vitesses du fluide sur les surfaces 
du contact s’êcrivent : 



fz = 0 u = Uj ; v = Vj et w = 0 

lz = h u = U 2 ; v = V 2 et w = W 2 




Fig. 7 : Système d'axes en coordonnêes cylindriques 
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